Diversity of major histocompatibility complex (MHC) and natural killer cell receptor (NKR) genes and their interactions in domestic horses.

The immunogenome is the part of the genome that underlies immune mechanisms and evolves under various selective pressures. Two complex regions of the immunogenome, major histocompatibility complex (MHC) and natural killer cell receptor (NKR) genes, play an important role in the response to selective pressures of pathogens. Their importance is expressed by their genetic polymorphism at the molecular level, and their diversity associated with different types of diseases at the population level. Findings of associations between specific combinations of MHC/NKR haplotypes with different diseases in model species suggest that these gene complexes did not evolve independently. No such associations have been described in horses so far. The aim of the study was to detect associations between MHC and NKR gene/microsatellite haplotypes in three horse breed groups (Camargue, African, and Romanian) by statistical methods; chi-square test, Fisher's exact test, Pearson's goodness-of-fit test and logistic regression. Associations were detected for both MHC/NKR genes and microsatellites; the most significant associations were found between the most variable KLRA3 gene and the EQCA-1 or EQCA-2 genes. This finding supports the assumption that the KLRA3 is an important receptor for MHC I and that interactions of these molecules play important roles in the horse immunity and reproduction. Despite some limitations of the study such as low numbers of horses or lack of knowledge of the selected genes functions, the results were consistent across different statistical methods and remained significant even after overconservative Bonferroni corrections. We therefore consider them biologically plausible.

Newly identified variability of the antigen binding site coding sequences of the equine major histocompatibility complex class I and class II genes.

The major histocompatibility complex (MHC) with its class I and II genes plays a crucial role in the immune response to pathogens by presenting oligopeptide antigens to various immune response effector cells. In order to counteract the vast variability of infectious agents, MHC class I and II genes usually retain high levels of SNPs mainly concentrated in the exons encoding the antigen binding sites. The aim of the study was to reveal new variability of selected MHC genes with a special focus on MHC class I physical haplotypes. Long-range NGS to was used to identify exon 2-exon 3 alleles in three genetically distinct horse breeds. A total of 116 allelic variants were found in the MHC class I genes Eqca-1, Eqca-2, Eqca-7 and Eqca-Ψ, 112 of which were novel. The MHC class II DRA locus was confirmed to comprise five exon 2 alleles, and no new sequences were observed. Additional variability in terms of 15 novel exon 2 alleles was identified in the DQA1 locus. Extensive overall variability across the entire MHC region was confirmed by an analysis of MHC-linked microsatellite loci. Both diversifying and purifying selection were detected within the MHC class I and II loci analyzed.

Microsatellite markers of the major histocompatibility complex genomic region of domestic camels.

We identified and characterized 11 polymorphic microsatellite markers suitable for routine testing (three in the MHC class I sub-region, four in MHC class II and four in the MHC class III sub-region) of dromedaries and Bactrian camels. In total, 38 dromedaries and 33 Bactrian camels were genotyped, and interspecific differences were observed in the numbers of alleles and in allelic frequencies, as well as in the observed heterozygosity. These loci may be used as markers to study the adaptive genetic diversity of the MHC region in Old World camels.

The natural cytotoxicity receptor genes in the family Felidae.

Natural killer (NK) cells belong to the innate immune system. The germline-encoded natural killer cell receptors represent activating and inhibitory receptors regulating multiple NK cell activities. The natural cytotoxicity receptors (NCRs) are activating natural cytotoxicity triggering receptors 1, 2, and 3 (NKp46, NKp44, and NKp30), encoded by the genes NCR1, NCR2, and NCR3, respectively. NCRs may be expressed in different cell types engaged in mechanisms of innate and adaptive immunity. The family Felidae, comprising the domestic cat and a wide variety of free-ranging species represents a well-suited model for biomedical and evolutionary studies. We characterized the NCR1, NCR2, and NCR3 genes in a panel of felid species. We confirmed the presence of potentially functional genes NCR1, NCR2, and NCR3 in all species. All three genes are conserved within the family and are similar to other phylogenetically related mammalian families. The NCR1 and NCR2 phylogenetic trees based on both nucleotide and protein sequences corresponded to the current zoological taxonomy, with some exceptions suggesting effects of different selection pressures in some species. Highly conserved NCR3 sequences did not allow a robust phylogenetic analysis. Most interspecific differences both at the nucleotide and protein level were found in NCR2. Within species, the most polymorphic CDS was detected in NCR1. Selection analyses indicated the effects of purifying selection on individual amino acid sites in all three genes. In stray cats, a rather high intraspecific diversity was observed.

Genomic analyses show extremely perilous conservation status of African and Asiatic cheetahs (Acinonyx jubatus).

We live in a world characterized by biodiversity loss and global environmental change. The extinction of large carnivores can have ramifying effects on ecosystems like an uncontrolled increase in wild herbivores, which in turn can have knock-on impacts on vegetation regeneration and communities. Cheetahs (Acinonyx jubatus) serve important ecosystem functions as apex predators; yet, they are quickly heading towards an uncertain future. Threatened by habitat loss, human-wildlife conflict and illegal trafficking, there are only approximately 7100 individuals remaining in nature. We present the most comprehensive genome-wide analysis of cheetah phylogeography and conservation genomics to date, assembling samples from nearly the entire current and past species' range. We show that their phylogeography is more complex than previously thought, and that East African cheetahs (A. j. raineyi) are genetically distinct from Southern African individuals (A. j. jubatus), warranting their recognition as a distinct subspecies. We found strong genetic differentiation between all classically recognized subspecies, thus refuting earlier findings that cheetahs show only little differentiation. The strongest differentiation was observed between the Asiatic and all the African subspecies. We detected high inbreeding in the Critically Endangered Iranian (A. j. venaticus) and North-western (A. j. hecki) subspecies, and show that overall cheetahs, along with snow leopards, have the lowest genome-wide heterozygosity of all the big cats. This further emphasizes the cheetah's perilous conservation status. Our results provide novel and important information on cheetah phylogeography that can support evidence-based conservation policy decisions to help protect this species. This is especially relevant in light of ongoing and proposed translocations across subspecies boundaries, and the increasing threats of illegal trafficking.

Comparative Genomics of the Major Histocompatibility Complex (MHC) of Felids.

This review summarizes the current knowledge on the major histocompatibility complex (MHC) of the family Felidae. This family comprises an important domestic species, the cat, as well as a variety of free-living felids, including several endangered species. As such, the Felidae have the potential to be an informative model for studying different aspects of the biological functions of MHC genes, such as their role in disease mechanisms and adaptation to different environments, as well as the importance of genetic diversity for conservation issues in free-ranging or captive populations. Despite this potential, the current knowledge on the MHC in the family as a whole is fragmentary and based mostly on studies of the domestic cat and selected species of big cats. The overall structure of the domestic cat MHC is similar to other mammalian MHCs following the general scheme "centromere-MHC class I-MHC class III-MHC class II" with some differences in the gene contents. An unambiguously defined orthologue of the non-classical class I HLA-E gene has not been identified so far and the class II DQ and DP genes are missing or pseudogenized, respectively. A comparison with available genomes of other felids showed a generally high level of structural and sequence conservation of the MHC region. Very little and fragmentary information on in vitro and/or in vivo biological functions of felid MHC genes is available. So far, no association studies have indicated effects of MHC genetic diversity on a particular disease. No information is available on the role of MHC class I molecules in interactions with Natural Killer (NK) cell receptors or on the putative evolutionary interactions (co-evolution) of the underlying genes. A comparison of complex genomic regions encoding NK cell receptors (the Leukocyte Receptor Complex, LRC and the Natural Killer Cell Complex, NKC) in the available felid genomes showed a higher variability in the NKC compared to the LRC and the MHC regions. Studies of the genetic diversity of domestic cat populations and/or specific breeds have focused mainly on DRB genes. Not surprisingly, higher levels of MHC diversity were observed in stray cats compared to pure breeds, as evaluated by DRB sequencing as well as by MHC-linked microsatellite typing. Immunogenetic analysis in wild felids has only been performed on MHC class I and II loci in tigers, Namibian leopards and cheetahs. This information is important as part of current conservation tasks to assess the adaptive potential of endangered wild species at the human-wildlife interface, which will be essential for preserving biodiversity in a functional ecosystem.

The Population Diversity of Candidate Genes for Resistance/Susceptibility to Coronavirus Infection in Domestic Cats: An Inter-Breed Comparison.

Feline coronavirus (FCoV) is a complex pathogen causing feline infectious peritonitis (FIP). Host genetics represents a factor contributing to the pathogenesis of the disease. Differential susceptibility of various breeds to FIP was reported with controversial results. The objective of this study was to compare the genetic diversity of different breeds on a panel of candidate genes potentially affecting FCoV infection. One hundred thirteen cats of six breeds were genotyped on a panel of sixteen candidate genes. SNP allelic/haplotype frequencies were calculated; pairwise FST and molecular variance analyses were performed. Principal coordinate (PCoA) and STRUCTURE analyses were used to infer population structure. Interbreed differences in allele frequencies were observed. PCoA analysis performed for all genes of the panel indicated no population substructure. In contrast to the full marker set, PCoA of SNP markers associated with FCoV shedding (NCR1 and SLX4IP) showed three clusters containing only alleles associated with susceptibility to FCoV shedding, homozygotes and heterozygotes for the susceptibility alleles, and all three genotypes, respectively. Each cluster contained cats of multiple breeds. Three clusters of haplotypes were identified by PCoA, two clusters by STRUCTURE. Haplotypes of a single gene (SNX5) differed significantly between the PCoA clusters.

Innate and Adaptive Immune Genes Associated with MERS-CoV Infection in Dromedaries.

The recent SARS-CoV-2 pandemic has refocused attention to the betacoronaviruses, only eight years after the emergence of another zoonotic betacoronavirus, the Middle East respiratory syndrome coronavirus (MERS-CoV). While the wild source of SARS-CoV-2 may be disputed, for MERS-CoV, dromedaries are considered as source of zoonotic human infections. Testing 100 immune-response genes in 121 dromedaries from United Arab Emirates (UAE) for potential association with present MERS-CoV infection, we identified candidate genes with important functions in the adaptive, MHC-class I (HLA-A-24-like) and II (HLA-DPB1-like), and innate immune response (PTPN4, MAGOHB), and in cilia coating the respiratory tract (DNAH7). Some of these genes previously have been associated with viral replication in SARS-CoV-1/-2 in humans, others have an important role in the movement of bronchial cilia. These results suggest similar host genetic pathways associated with these betacoronaviruses, although further work is required to better understand the MERS-CoV disease dynamics in both dromedaries and humans.

Crimean-Congo Hemorrhagic Fever Virus Past Infections Are Associated with Two Innate Immune Response Candidate Genes in Dromedaries.

Dromedaries are an important livestock, used as beasts of burden and for meat and milk production. However, they can act as an intermediate source or vector for transmitting zoonotic viruses to humans, such as the Middle East respiratory syndrome coronavirus (MERS-CoV) or Crimean-Congo hemorrhagic fever virus (CCHFV). After several outbreaks of CCHFV in the Arabian Peninsula, recent studies have demonstrated that CCHFV is endemic in dromedaries and camel ticks in the United Arab Emirates (UAE). There is no apparent disease in dromedaries after the bite of infected ticks; in contrast, fever, myalgia, lymphadenopathy, and petechial hemorrhaging are common symptoms in humans, with a case fatality ratio of up to 40%. We used the in-solution hybridization capture of 100 annotated immune genes to genotype 121 dromedaries from the UAE tested for seropositivity to CCHFV. Through univariate linear regression analysis, we identified two candidate genes belonging to the innate immune system: FCAR and CLEC2B. These genes have important functions in the host defense against viral infections and in stimulating natural killer cells, respectively. This study opens doors for future research into immune defense mechanisms in an enzootic host against an important zoonotic disease.

Candidate Gene Markers Associated with Fecal Shedding of the Feline Enteric Coronavirus (FECV).

The Feline coronavirus (FCoV) can cause a fatal disease, the Feline Infectious Peritonitis. Persistent shedders represent the most important source of infection. The role of the host in FCoV fecal shedding is unknown. The objective of this study was to develop gene markers and to test their associations with FCoV shedding patterns. Fecal samples were taken from 57 cats of 12 breeds on the day 0 and after 2, 4 and 12 months. Variation from persistent and/or high-intensity shedding to no shedding was observed. Thirteen immunity-related genes were selected as functional and positional/functional candidates. Positional candidates were selected in a candidate region detected by a GWAS analysis. Tens to hundreds of single nucleotide polymorphisms (SNPs) per gene were identified using next generation sequencing. Associations with different phenotypes were assessed by chi-square and Fisher's exact tests. SNPs of one functional and one positional candidate (NCR1 and SLX4IP, respectively) and haplotypes of four genes (SNX5, NCR2, SLX4IP, NCR1) were associated with FCoV shedding at pcorected < 0.01. Highly significant associations were observed for extreme phenotypes (persistent/high-intensity shedders and non-shedders) suggesting that there are two major phenotypes associated with different genotypes, highly susceptible cats permanently shedding high amounts of viral particles and resistant non-shedders.

The Camel Adaptive Immune Receptors Repertoire as a Singular Example of Structural and Functional Genomics.

The adaptive immune receptors repertoire is highly plastic, with its ability to produce antigen-binding molecules and select those with high affinity for their antigen. Species have developed diverse genetic and structural strategies to create their respective repertoires required for their survival in the different environments. Camelids, until now, considered as a case of evolutionary innovation because of their only heavy-chain antibodies, represent a new mammalian model particularly useful for understanding the role of diversity in the immune system function. Here, we review the structural and functional characteristics and the current status of the genomic organization of camel immunoglobulins (IG) or antibodies, α/ß and γ/δ T cell receptors (TR), and major histocompatibility complex (MHC). In camelid humoral response, in addition to the conventional antibodies, there are IG with "only-heavy-chain" (no light chain, and two identical heavy gamma chains lacking CH1 and with a VH domain designated as VHH). The unique features of these VHH offer advantages in biotechnology and for clinical applications. The TRG and TRD rearranged variable domains of Camelus dromedarius (Arabian camel) display somatic hypermutation (SHM), increasing the intrinsic structural stability in the γ/δ heterodimer and influencing the affinity maturation to a given antigen similar to immunoglobulin genes. The SHM increases the dromedary γ/δ repertoire diversity. In Camelus genus, the general structural organization of the TRB locus is similar to that of the other artiodactyl species, with a pool of TRBV genes positioned at the 5' end of three in tandem D-J-C clusters, followed by a single TRBV gene with an inverted transcriptional orientation located at the 3' end. At the difference of TRG and TRD, the diversity of the TRB variable domains is not shaped by SHM and depends from the classical combinatorial and junctional diversity. The MHC locus is located on chromosome 20 in Camelus dromedarius. Cytogenetic and comparative whole genome analyses revealed the order of the three major regions "Centromere-ClassII-ClassIII-ClassI". Unexpectedly low extent of polymorphisms and haplotypes was observed in all Old World camels despite different geographic origins.

The Major Histocompatibility Complex of Old World Camels-A Synopsis.

This study brings new information on major histocompatibility complex (MHC) class III sub-region genes in Old World camels and integrates current knowledge of the MHC region into a comprehensive overview for Old World camels. Out of the MHC class III genes characterized, TNFA and the LY6 gene family showed high levels of conservation, characteristic for MHC class III loci in general. For comparison, an MHC class II gene TAP1, not coding for antigen presenting molecules but functionally related to MHC antigen presenting functions was studied. TAP1 had many SNPs, even higher than the MHC class I and II genes encoding antigen presenting molecules. Based on this knowledge and using new camel genomic resources, we constructed an improved genomic map of the entire MHC region of Old World camels. The MHC class III sub-region shows a standard organization similar to that of pig or cattle. The overall genomic structure of the camel MHC is more similar to pig MHC than to cattle MHC. This conclusion is supported by differences in the organization of the MHC class II sub-region, absence of functional DY genes, different organization of MIC genes in the MHC class I sub-region, and generally closer evolutionary relationships of camel and porcine MHC gene sequences analyzed so far.

The major histocompatibility complex of Old World camelids: Class I and class I-related genes.

The genomic structure of the Major Histocompatibility Complex (MHC) region and variation in selected MHC class I related genes in Old World camels, Camelus bactrianus and Camelus dromedaries were studied. The overall genomic organization of the camel MHC region follows a general pattern observed in other mammalian species and individual MHC loci appear to be well conserved. Selected MHC class I genes B-67 and BL3-7 exhibited unexpectedly low variability, even when compared to other camel MHC class I related genes MR1 and MICA. Interspecific SNP and allele sharing are relatively common, and frequencies of heterozygotes are usually low. Such a low variation in a genomic region generally considered as one of the most polymorphic in vertebrate genomes is unusual. Evolutionary relationships between MHC class I related genes and their counterparts from other species seem to be rather complex. Often, they do not follow the general evolutionary history of the species concerned. Close evolutionary relationships of individual MHC class I loci between camels, humans and dogs were observed. Based on the results of this study and on our data on MHC class II genes, the extent and the pattern of polymorphism of the MHC region of Old World camelids differed from most mammalian groups studied so far. Camels thus seem to be an important model for our understanding of the role of genetic diversity in immune functions, especially in the context of unique features of their immunoglobulin and T-cell receptor genes.

The major histocompatibility complex in Old World camelids and low polymorphism of its class II genes.

BACKGROUND: The Major Histocompatibility Complex (MHC) is a genomic region containing genes with crucial roles in immune responses. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. To counteract the high variability of pathogens, the MHC evolved into a region of considerable heterogeneity in its organization, number and extent of polymorphism. Studies of MHCs in different model species contribute to our understanding of mechanisms of immunity, diseases and their evolution. Camels are economically important domestic animals and interesting biomodels. Three species of Old World camels have been recognized: the dromedary (Camelus dromedarius), Bactrian camel (Camelus bactrianus) and the wild camel (Camelus ferus). Despite their importance, little is known about the MHC genomic region, its organization and diversity in camels. The objectives of this study were to identify, map and characterize the MHC region of Old World camelids, with special attention to genetic variation at selected class MHC II loci. RESULTS: Physical mapping located the MHC region to the chromosome 20 in Camelus dromedarius. Cytogenetic and comparative analyses of whole genome sequences showed that the order of the three major sub-regions is "Centromere - Class II - Class III - Class I". DRA, DRB, DQA and DQB exon 2 sequences encoding the antigen binding site of the corresponding class II antigen presenting molecules showed high degree of sequence similarity and extensive allele sharing across the three species. Unexpectedly low extent of polymorphism with low numbers of alleles and haplotypes was observed in all species, despite different geographic origins of the camels analyzed. The DRA locus was found to be polymorphic, with three alleles shared by all three species. DRA and DQA sequences retrieved from ancient DNA samples of Camelus dromedarius suggested that additional polymorphism might exist. CONCLUSIONS: This study provided evidence that camels possess an MHC comparable to other mammalian species in terms of its genomic localization, organization and sequence similarity. We described ancient variation at the DRA locus, monomorphic in most species. The extent of molecular diversity of MHC class II genes seems to be substantially lower in Old World camels than in other mammalian species.

Flow cytometry-based enumeration and functional characterization of CD8 T regulatory cells in patients with multiple myeloma before and after lenalidomide plus dexamethasone treatment.

BACKGROUND: Multiple myeloma (MM) is a malignancy of plasma cells frequently associated with immune abnormalities. Several studies have confirmed that in MM immune deregulation can be mediated by increased numbers of CD4 T regulatory (Treg) cells, and these cells were also associated with poor outcome. In this study, we aimed to study CD8 Treg cells before and after lenalidomide plus dexamethasone (len-dex) treatment in MM patients. METHODS: Using flow cytometry, we enumerated and assessed suppressive function of CD8 Treg cells in 16 MM patients before and after len-dex treatment. RESULTS: Numbers of CD8 Treg cells (CD8+CD25hi+FoxP3+) (P < 0.01) were significantly increased in MM patients (before treatment) compared to healthy donors. However, no significant changes were observed in CD4 and CD8 T cells. A significant increase in CD8 Treg cells was observed after len-dex treatment compared to pre-treatment but no significant difference was observed in CD4 and CD8 T cells. Proliferation assay data showed that CD8 Treg cells inhibited proliferation of CD4 T cells and IFN-γ secretion in a concentration dependent manner. Suppressive activity of CD8 Treg cells did not differ significantly between healthy donors, untreated and len-dex treated MM patients. A significant abnormal level of IL-10 was observed from proliferation assays of untreated and len-dex treated MM patients compared to healthy donors (P ≤ 0.03). CONCLUSIONS: Using flow cytometry, we have shown that suppressive CD8 Treg cells are increased in MM patients and len-dex treatment is unable to control these suppressive CD8 Treg cells.

Flow cytometry based enumeration and functional characterization of CD8 T regulatory cells in patients with multiple myeloma before and after lenalidomide plus dexamethasone treatment.

BACKGROUND: Multiple myeloma (MM) is a malignancy of plasma cells frequently associated with immune abnormalities. Several studies have confirmed that in MM immune deregulation can be mediated by increased numbers of CD4 T regulatory (Treg) cells, and these cells were also associated with poor outcome. In the present study, we aimed to study CD8 Treg cells before and after lenalidomide plus dexamethasone (len-dex) treatment in MM patients. METHODS: Using flow cytometry, we enumerated and assessed suppressive function of CD8 Treg cells in 16 MM patients before and after len-dex treatment. RESULTS: Numbers of CD8 Treg cells (CD8+CD25hi+FoxP3+) (P<0.01) were significantly increased in MM patients (before treatment) compared to healthy donors. However, no significant changes were observed in CD4 and CD8 T cells. A significant increase in CD8 Treg cells was observed after len-dex treatment compared to pre-treatment but no significant difference was observed in CD4 and CD8 T cells. Proliferation assay data showed that CD8 Treg cells inhibited proliferation of CD4 T cells and IFN-γ secretion in a concentration dependent manner. Suppressive activity of CD8 Treg cells did not differ significantly between healthy donors, untreated and len-dex treated MM patients. A significant abnormal level of IL-10 was observed from proliferation assays of untreated and len-dex treated MM patients compared to healthy donors (P≤0.03). CONCLUSION: Using flow cytometry, we have shown that suppressive CD8 Treg cells are increased in MM patients and len-dex treatment is unable to control these suppressive CD8 Treg cells. © 2013 Clinical Cytometry Society.